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摘要 : 叶 面 积 指 数 ( Leaf Area Index，LAI) 高度 综 合 了 植被 水 平 履 新 状况 和 垂直 结构 ;以 及 枯 校 落叶 层 厚 注 和 地 下 生物 量 多 少 ， 
是 植被 影响 土壤 侵蚀 的 主要 方面 。 区 域 尺度 的 时 间 序 列 叶 面 积 指 数 揭示 了 区 域 土 壤 侵蚀 的 演化 过 程 因此 ,及 时 准确 地 掌握 
区 域 尺度 上 长 时 间 序 列 的 植被 LAI, 对 人 研究 土壤 侵蚀 动态 变化 与 植被 的 关系 至 关 重 要 6 选择 南京 市 1988 一 2013 年 10 期 遥感 影 
像 ,基于 反问 传播 ( Back Propagation，BP ) 神 经 网 络 构建 LAI 反 演 模型 ,进行 了 长 时 间 序 列 的 叶 面 积 指数 反 演 。 结 合 2009 和 
2010 年 LAI 实测 值 ,验证 与 探讨 了 该 模型 的 评价 精度 与 适应 性 。 结 果 表 明 :;( 1) 该 模型 拟 合 度 较 高 ,2009 和 2010 年 平均 相对 误 
差 . 均 方 根 误差 .相关 系数 分 别 是 0.2395 和 0.2174, 0.2962 和 0.2581, 0:7713 和 '0.6844 ;各 项 精度 评价 指标 均 较 好 ; (2) 统 计 分 
析 去 除 耕 地 后 全 市 LAI 变化 , 低 植被 覆盖 (LAI<2) 面积 不 断 增 加 ,高 植被 覆盖 区 (LAI>3) 面 积 先 减少 后 增加 ,耕地 面积 不 断 减 
少 ,符合 南京 市 的 发 展 变化 规律 ;(3) 主 城区 LAI 年 际 变化 与 其 他 学 者 得 到 的 南京 市 植被 盖 度 变化 趋势 一 致 , 反 演 结果 的 时 序 
性 较 高 。 本 文 提出 的 基于 反问 传播 神经 网 络 模 型 反 演 长 时 间 序 列 LAT 是 可 行 的 ,为 区 域 尺度 土壤 侵蚀 定量 遥感 监测 提供 新 
途径 。 
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Abstract: Leaf Area Index( LAI) can reflect the horizontal coverage, vertical structure of vegetation, the thickness of the 
litter layer arid-the amount of underground biomass, which ls the main aspect of vegetation affecting soil erosion. lt is very 
important to monitor the changes in the amount of soil erosion, for useful information to guide the planning of soil and water 
conservati6n, protect the soil and water resources and control the soil erosion. Therefore, the method by which we obtain 
high quality and long sequential LAIL at a regional scale is very important for analyzing the relationship between the dynamic 
changes in soil erosion and vegetation. Previous studies showed that the neural network had an incomparable superiority in 
terms of complex, nonlinear data fitting and pattern recognition, and had been successfully applied to inverse the LAI in 


Nanjing based on the multi-spectral remote sensing data derived from the Landsat 8 Operational Land Imager( OLI), four 
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types of vegetation indices ( Normalized Difference Vegetation Index, NDVI; Ratio Vegetation Index, RVI; Soil-adjusted 
Vegetation Index, SAVI; Modified Soil adjusted Vegetation Index, MSAVI), and measured LAI data. The results showed 
that the accuracy of retrieval was good. In this paper, we used the Back Propagation( BP) neural network model to inverse 
the LAI in Nanjing during 1988 一 2013 based on the data derived from Landsat 8 OLI and Landsat 3 Thematic Mapper 
(TM ). Based on the measured values of LAI in 2009 and 2010, the evaluation accuracy and adaptability of the model were 
verified and discussed. The results showed that:(1) The model had a fitting of higher degree, and average relative errors 
(MAPE ) , root mean square errors( RMSE ) , and correlation coefficients( R) in 2009 and 2010 of 0.2395 and 0.2174, 
0.2962 and 0.2581, and 0.7713 and 0.6844, respectively. Each accuracy evaluation index was good. (2) Following 
statistical analysis, we found that the low vegetation coverage area( LAT < 2) exhibited an increasing trend, “the high 
vegetation coverage area( LAT > 3) presented a first decreasing and then increasing trend, while the cultivated Jand area 
decreased with the rapid development in Nanjing.(3) To analyze accurately the LAI, we extracted the LAI "in the main 
urban area, and found that there was a relatively high inversion value, and the inter-annual’change im LAITwas consistent 
with the change in vegetation coverage in Nanjing reported by other studies. Therefore ，we could see that the BP neural 
network model had a high accuracy for the time series LAIT inversion. It provides a new way for qtiantitative remote sensing 
monitoring of regional soil erosion. Moreover, because of other potential limiting factors ,such as the errors produced by the 
BP neural network model the large area of the inversion area, the complexity of vegetation types and community structure, 
etc., the inversion accuracy of LAIL through remote sensing still needs to be exploredy and the inversion method improved. 
We will try to establish a multi-angle LAI inversion method to corstruct the coupling model of LAI and soil erosion or 


quantitative fusion of multi-source remote sensing images in the future study: 


Key Words: Nanjing City; BP neural network; long time-series;y Leaf Area Index( LAI1) ; soil erosion 


叶 面 积 指数 (Leaf Area Index，LAI) 是 植被 冠 层 的 重要 结构 参数 ,是 陆 面 生态 过 程 模型 的 重要 输入 参量 ， 
对 植被 冠 层 截流 .光合 .呼吸 . 蒸 散 、 降 水 堆 持 、 碳 循环 和 土壤 侵蚀 等 起 到 关键 作用 , 它 能 反映 出 植被 的 水 平 履 
盖 状 次 和 往 直 结 构 ,甚至 枯 术 落叶 层 郁 薄 和 地 下 生物 量 多 少 ,这 些 属 性 正 是 植被 影响 土壤 侵蚀 的 主要 方 
面 … 。 杨 勤 科 等 认为 LAI 更 适合 代替 植被 覆盖 度 作 为 水 土 保持 定量 评价 的 重要 指标 …”。 王 库 等 利用 卫星 
影像 数据 对 江西 省 兴国 县 区 域 的 植被 叶 面 积 指 数 分 布 与 土壤 侵蚀 分 布 关系 进行 了 分 析 , 定 性 地 说 明 植 被 叶 面 
积 指数 与 土壤 侵蚀 有 着 良好 的 相关 性 "”。 孙 佳 佳 等 研究 了 马尾 松 纯 林 马尾 松林 草 的 不 同 覆盖 度 和 植被 叶 
面积 指数 与 土壤 侵蚀 的 关系 ,结果 表明 选择 LAI 表征 土壤 侵蚀 模 数 来 评价 植被 水 土 保持 效益 ,更 为 稳定 和 可 
靠 ” 。 林 态 等 基于 工 AI 构 建 了 USLE 方程 中 的 植被 覆盖 管理 措施 C 因子 ,结果 显示 LAI 取代 植被 覆盖 度 进行 
水 土 保持 定量 估算 的 精度 得 到 明显 提高 ”。 因 此 ,区 域 尺度 长 时 间 序 列 LAI 遥感 反 演 的 精度 直接 影响 土壤 
侵蚀 定量 评价 。 

遥感 技术 具有 大 范围 .快速 和 连续 时 间 序 列 等 优点 ,是 估算 区 域 乃 至 全 球 太 度 时 间 序 列 LAI 的 重要 技术 
手段 ,其 方法 有 统计 模型 法 .物理 模型 方法 等 多 种 ”。 统 计 模型 法 简单 快捷 ,被 很 多 学 者 应 用 于 LAI 反 演 或 
LAT 反 演 方法 对 比 研 究 中 ,但 由 于 统计 模型 易 受 土壤 背景 等 多 种 外 在 因素 影响 , 反 演 精度 往往 不 高 ,缺少 可 移 
植 性 ”。 物 理 模 型 方法 的 理论 基础 完善 ,模型 的 参数 具有 明确 的 物理 意义 ,并 可 对 作用 机 理 进行 适当 的 数 
学 描述 ,但 此 类 模型 一 般 是 非 线 性 的 ,输入 参数 多 ,方程 复杂 ,适用 性 较 差 ,日 对 非 主 要 因素 有 过 多 的 忽略 或 假 
定 “"”。 神 经 网 络 在 对 复杂 、 非 线性 数据 拟 合 及 模式 识别 方面 有 着 无 可 比拟 的 优势 “"” , 杨 敏 等 基于 2013 
年 Landsat 8 Operational Land Imager( Landsat 8 OLI) 数据 构建 了 基于 后 向 反射 (Back Propagation, BP ) 神经 网 
络 的 LAI 反 演 模 型 , 反 演 精度 明显 高 于 非 线性 回归 模型 ,模型 具有 较 高 的 空间 可 靠 性 。 

近年 来 , 叶 面 积 指数 反 演 在 农业 方面 已 做 了 大 量 研 究 ”” ,但 是 在 林业 方面 研究 相对 较 少 。 搭 载 在 Terra 
和 Aqua 两 条 卫星 上 的 Moderate-resolution Imaging Spectroradiometer( MODIS ) 传 感 融 能 够 长 时 间 收 集 全 球 陆地 
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表面 LAI 的 变化 信息 。 然 而 ,目前 发 布 的 MODIS LAI 数据 产品 的 时 空 不 连续 性 和 较 低 的 空间 分 状 率 制约 着 
MODIS LAI 产品 在 土壤 侵蚀 模型 方面 的 应 用 ,高 分 辨 率 长 时 间 序 列 的 LAI 应 用 到 土壤 侵蚀 方面 研究 相对 较 
少 。 因 此 本 文 以 南京 市 为 研究 区 域 ,应 用 BP 神经 网 络 模型 ,基于 1988 一 2013 年 Landsat TM 遥感 数据 , 反 演 
长 时 间 序 列 的 叶 面 积 指数 ,对 比分 析 LAI 的 年 际 变化 与 植被 覆 善 度 变化 的 关系 ,为 区 域 土 壤 侵蚀 动态 变化 遥 
感 监 测 提供 新 方法 。 


1 研究 区 概况 与 研究 方法 


1.1 研究 区 概况 

南京 市 位 于 长 江 下 游 中 部 地 区 ,江苏 省 西南 部 , 介 于 31°14’ 一 32°37’'N,118?22' 一 119°14'E 之 间 。 全 市 行 
政 区 域 总 面积 6587.02 km ,地形 以 低 山 丘陵 为 骨架 ,以 环 状 山 、 条 带 山 . 自 状 盆地 为 主要 特色 。 和 气候 属于 北 
亚热带 季风 气候 , 四季 分 明 ,雨量 充沛 ,年 平均 气温 16% ,年 平均 降雨 117 d ,年 平均 降雨 量 1106 mm ,无 霜 期 
237 d。 南 京 地 区 的 土壤 在 北 、 中 部 广大 地 区 为 黄 标 壤 ( 地 市 性 土壤 ) ,南部 与 安徽 省 接壤 人 处 有 小 面积 的 红壤 。 
植被 类 型 属 常 绿 落叶 阔 叶 混 交 林 。 南 京 地 区 人 口 密集 ,属于 农业 活动 强烈 区 ,自然 植被 在 历史 上 屡 遭 严重 破 
坏 ,几乎 全 部 消失 , 现 有 植被 多 属 次 生性 质 ,其 中 人 工 林 面积 大 于 自然 恢复 的 次 生 林 。 境 内 现 有 林业 用 地 约 
840 km?。 用 材 林 和 生态 林 约 570 km2 ,经 济 林 和 人 竹林 210 km2z 。 用 材 林 和 生态 林 的 树种 主要 有 马尾 松 、 黑 松 、 
杉木 .国外 松 、 麻 栎 .刺槐 .水池 杉 、 柏 类 等 。 经 济 林 以 茶 果 和 桑 为 主 。 人 竹林 以 毛竹 为 主 ,集中 于 丘陵 山区 。 
1.2 ”遥感 影像 数据 及 处 理 


1) 遥感 影像 数据 
本 文 用 于 反 演 的 遥感 影像 数据 为 Landsat 8 0ET 及 Eandsat S$_Thematic Mapper( Landsat 5 TM ) 数据 , 共 10 


景 影像 数据 。 影 像 选取 要 求 为 :选取 时 间 为 当年 7 一 9 月份 ,是 植被 生长 最 为 成 盛 的 月 份 ,也 是 降雨 最 为 集中 ， 
最 容易 发 生 水 土 流失 的 月 份 , 为 了 保证 反 演 精度 ?研究 区 还 要 人 太 量 选取 无 云 履 盖 的 高 质量 影像 ,但 是 2000、 
2004 及 2009 年 在 7 一 9 月 没有 符合 质量 要 求 的 影像 , 故 选 择 了 10 月 的 影像 。 最 终 影像 数据 如 表 1 所 示 。 

Landsat 5 为 光学 对 地 观测 卫星 ,是 美国 陆地 卫星 系列 ( Landsat 卫星 ) 的 第 五 颗 卫 星 , 于 1984 年 3 月 1 日 
发 射 ,2011 年 11 月 18 日 ,美国 地 质 调查 局 发 布 消息 称 由 于 卫星 上 的 放大 带 迅 速 老 化 ,日 前 已 停止 获取 
Landsat 5 的 卫星 遥感 影像 。Landsat8_ OLIT 多 光谱 遥感 影像 数据 由 美国 地 质 勘 探 局 提供 。OLI 包括 了 ETM+ 传 
感 带 所 有 的 波段 ,为 了 避免 大 气 吸收 特征 , OLI 对 波段 进行 了 重新 调整 ,其 中 比较 大 的 调整 是 OLI Band5 
(0.845 一 0.885 pm) ,排除 了 0.825 pm 处 水 汽 吸 收 特 征 。 


表 工 用 于 南京 市 时 间 序 列 LAI 反 演 的 遥感 影像 数据 统计 


Table 1 The remote sensing image data used to time series LAI inversion of Nanjing City 


序号 Number 成 像 时 间 Imaging time 传感器 Sensor 分 辨 率 Resolution/m 轨道 号 Path/row 
1 1988-07-05 LANDSAT 5_TM 30 120/38 
2 1994-07-22 LANDSAT 5_TM 30 120/38 
3 2000- 10- 10 LANDSAT 5_TM 30 120/38 
4 2002-07- 12 LANDSAT 5_TM 30 120/38 
5 2004- 10-21 LANDSAT 5_TM 30 120/38 
6 2006-08-08 LANDSAT 5_TM 30 120/38 
7 2007-07-26 LANDSAT 5_TM 30 120/38 
8 2009- 10-03 LANDSAT 5_TM 30 120/38 
9 2010-08- 19 LANDSAT 5_TM 30 120/38 
10 2013-08-11 LANDSAT_8_OLI 30 120/38 


2) 其 他 数据 源 包 括 ; 南 京 市 1986 .1996 .2002 年 三 期 土地 利用 / 禾 被 数据 ( 中科院 提供 ) ,2013 年 土地 利用 
窗 盖 数据 是 以 2013 年 Landsat8 0LI 数据 为 基础 ,参考 GOOGLE EARTH 高 分 辩 率 影像 和 2002 年 的 土地 利用 
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覆 被 数据 以 及 国家 科技 基础 条 件 平台 一 一 地 球 系统 科学 数据 共享 平台 获取 的 2013 年 土地 利用 覆 被 数据 进行 
人 机 交互 解 译 获得 ,经 实地 采样 验证 , 解 译 精度 在 88% 以 上 。 

3) 数 据 预 处 理 

大 气 校正 是 遥感 影像 获取 地 面 真实 反射 率 必 不 可 
少 的 步骤, 对 植被 定量 遥感 尤为 重要 “” 。 本 文采 用 
Fast line-of-sight atmospheric analysis of spectral 
hypercubes(FLAASH ) 模型 进行 大 气 校正 。 同 时 ,以 校 
正好 的 南京 市 2007 年 TM 遥感 影像 为 标准 底 图 ,采用 
二 次 多 项 式 拟 合法 ,对 影像 进行 几何 精 纠 正 ,误差 控制 
在 0.5 个 像 元 内 。 最 后 ,用 南京 边界 矢量 数据 进行 裁 
切 ,得 到 研究 区 域 10 景 影像 数据 。 预 处 理 流程 图 如 
ll 
1.3 LAI 野外 实测 数据 

为 了 验证 神经 网 络 反 演 模型 对 于 时 间 序 列 LAI 反 
演 精度 的 有 效 性 ,本 文 利用 2010 年 和 2009 年 的 LAI 实 
测 值 进行 反 演 模型 的 精度 验证 。 

LAI 数据 的 测定 是 采用 美国 LI-COR LAI- 2200 植 
物 冠 层 分 析 仪 。2009 年 , 样 点 数 19 个 , 主要 分 布 在 铜 
山林 场 。2010 年 , 样 点 数 37 个 , 主要 分 布 在 江宁 区 和 图 1 多 光谱 遥感 影像 预 处 理 流程 图 
铜 山林 场 。 所 有 测定 时 间 都 选择 在 6:30 一 9:00 点 之 间 Fig.1 Miilti-spectral remote sensing image pre-processing flow 
或 16:30 一 19:00 之 间 , 尽 量 避 免 因 太阳 光线 的 下 射 而 、 shart 
引起 的 测试 误差 ,每 个 样 地 分 别 在 4 个 角 和 中 心 位 置 各 
测量 一 次 , 取 5 次 均值 作为 结果 。 鉴 于 TM 影像 30m 的 空间 分 辩 率 ,采样 间距 均 大 于 30 mm, 每 个 样 点 均 由 
GARMIN 手持 GPS 接收 机 定位 ,坐标 系 为 WGS-84, 各 样 点 重复 测量 2 次 , 取 均 值 作为 结 
1.4 LAI 反 演 算法 

Landsat 8 OLI 多 光谱 遥感 影像 数据 包括 了 Enhanced Thematic Mapper(ETM+ ) 传 感 硕 所 有 的 波段 ,只 是 对 
部 分 波段 进行 了 调整 ,增加 了 两 个 波段 : 蓝 波 段 (Bandl ;0.433 一 0.453 ) 和 短波 红外 (Band9;1.360 一 1.390) ,其 
中 比较 大 的 调整 是 OLL Band5(0.845 一 0.885 pm) ,排除 了 0.825 pm 处 水 汽 吸 收 特征 。 因 此 ,基于 2013 年 OLI 
数据 所 建 的 模型 适用 于 前 25 年 的 TM-LAI 估算 。 本 文 应 用 隐 含 层 为 2 层 的 BP 神经 网 络 模型 对 1988 一 2013 
年 南京 市 LAT 和 进行 反 演 5 ,利用 2009 年 和 2010 年 LAI 的 实测 数据 检验 反 演 精度 。 但 是 由 于 传感器 的 差异 、 
大 气 等 因素 的 影响 ,每 景 影像 相同 地 物 的 反射 率 会 有 一 定 的 差异 ,因此 本 文 将 神经 网 络 的 训练 和 模拟 的 归 一 
化 部 分 做 了 相应 调整 , 即 :对 每 景 分 别 做 归 一 化 处 理 , 减 小 反射 率 差异 对 反 演 精度 造成 的 影响 。 
1S$， 精 度 评 价 方法 

平均 相对 误差 ( MAPE) 、 均 方 根 误差 (RMSE) 和 相关 系数 (R) 被 用 来 衡量 和 刻画 模型 的 精度 。 其 计算 公 
式 分 别 为 (公式 1 一 3): 


辐射 定 标 


几何 精 校正 


2007 年 TM 基准 影像 


南京 市 边界 矢量 数据 


> NiCBP) Xi 
i=1 Xi 
MAPE = 本 (1) 
> (wenes 一 Xi ) 
RMSE = | 一 6 


Nn 
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> (Wicppy — %) (Xi 一 无) 
CO (3) 


之 Cos -1 /> (a 2 
式 中 , xicap) 为 模型 反 演 值 LAI( BP);%; 为 LAI 实测 值 ; x 为 LAI 实测 值 的 均值 ; o” 为 样本 方差 ;n 为 样本 
个 数 。 
2 结果 与 分 析 


2.1 基于 地 面 实测 LAI 数据 的 误差 分 析 
为 了 进一步 验证 时 间 序 列 LAI 反 演 的 精度 ,本 文 分 别 利 用 2009 年 和 2010 年 的 LAI 实测 值 进行 检验 , 拟 
合 结果 如 图 2 ,模拟 精度 如 表 2。 
国 LAI 实 测 值 


2009 
4.0 国 LAI 反 演 值 


3.5 


3.0 


一 
Cn 
[LI | | 


叶 面 积 指数 LAI Leaf area index 
[>] 
Le 
S 
好 


3 

| 
| 
| 


234 567 891011 1213141516171819 2021 22232425262728 293031323334353637 
样本 点 号 The sample dot 


图 2 2009 年 和 2010 年 LAI 反 演 值 和 实测 值 
Fig.2 Fitting results of measured LAI and retrieval LAI in 2009 and 2010 


表 2 2009 年 和 2010 年 LAI 反 演 精度 
Table 2 LAI retrieval accuracy in 2009 and 2010 


采样 时 间 样本 点 数 


MAPE RMSE R 
Sampling time The sample dot 
2009 19 0.2395 0.2962 0.7713 
2010 37 0.2174 0.2581 0.6844 
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图 2 显示 模型 能 较 好 地 模拟 2009 年 和 2010 年 的 LAI, 模 拟 值 和 实测 值 吻 合 度 较 高 。 由 表 2 可 知 ,2009 
和 2010 年 平均 相对 误差 (MAPE ) . 均 方 根 误差 (RMSE) .相关 系数 (R) 分 别 是 0.2395 和 0.2174 .0.2962 和 
0.2581 .0.7713 和 0.6844 ,各 项 精度 评价 指标 均 较 好 。 

利用 基于 2013 年 Landsat8 OIL 数据 训练 生成 的 BP 神经 网 络 模型 ,成 功 模拟 了 1988 一 2013 年 的 LAI, 生 
成 1988 一 2013 年 LAI 分 布 图 ,如 图 3。 


1988 年 7 月 5 号 1994 年 7 月 22 号 2000 年 10 月 10 号 2002 年 7 月 12 号 


OO 


3 1988 一 2013 年 南京 市 LAI 分 布 图 
Fig.3 LAI distribution map of Nanjing City in 1988 一 2013 


从 图 3 中 可 以 直观 地 看 出 ,LAI 小 于 2 的 主要 分 布 在 建筑 物 比较 密集 .植被 稀 焉 的 城区 ,包括 江宁 区 北 
部 六合 区 南部 浦口 区 东南 部 栖霞 区 以 及 市 辖区 ;LAI 在 2 一 3 之 间 的 ,主要 分 布 在 建筑 物 较 多 、 有 部 分 植被 
的 城区 与 乡镇 ,面积 较 大 ,分 布 比较 分 散 ;LAI 在 3 一 4 之 间 的 ,主要 分 布 在 建筑 物 较 少 、 植 被 较 多 的 郊区 以 及 
农田 ,面积 最 大 ;LAI 大 于 4 的 主要 分 布 在 紫金 山 、 老 山 、 磊 府 山 将 军 山 、 铜 山 等 植被 比较 密集 的 山区 。 模 拟 
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结果 与 南京 市 的 植被 分 布 实际 情况 相符 。 
2.2 ”南京 市 植被 变化 分 析 

由 于 南京 市 四 周 有 大 量 农田 分 布 ,农田 受 人 为 影响 较为 明显 ,LAI 的 变化 不 仅 与 农田 的 面积 变化 有 关 , 农 
作物 的 长 势 也 将 直接 影响 LAI 的 大 小 ,为 了 排除 农作物 对 LAI 变化 的 影响 , 本文 根据 南京 市 1988、1994、 
2002 .2013 四 年 的 土地 利用 分 类 ,用 耕地 进行 掩 膜 处 理 , 生 成 去 除 耕 地 后 南京 市 的 LAI 分布 图 (图 4)。 


1988 年 


1994 年 2002 年 2013 年 


0 20km 
图 例 
各 区 县 边界 <2.0 mm 2.5 一 3.0 mm 3.5 一 4.0 
南京 水 域 国 2.0 一 2.5 mm 3.0 一 3.5 mm >4.0 


图 4 去 除 耕 地 后 南京 市 LAI 分 布 图 
Fig.4 LATI distribution map of Nanjing City without cultivated land 


如 图 4 所 示 , 可 以 看 出 南京 近 40 年 来 ,去 除 耕 地 后 建筑 用 地 等 面积 逐年 增加 ,由 主 城区 向 四 周 扩 散 , 相 应 


的 耕地 面积 逐渐 较 小 。 为 了 定量 地 了 解 LAL 的 变化 趋势 ,本 研究 统计 分 析 了 去 除 耕 地 后 全 市 的 低 LAI 区 域 
(LAI<2) 和 高 LAI 区 域 (LAI>3) 的 变化 趋 热 ( 图 5 六 


900 60 - 
0 0 3x sof 
外 600 -e-LAI>3 六 0 40 I 
pS 500 吕 二 这 
< 400 喘 8330r 
习 300 _。 拖 昌 < 20F 
鸯 ”200 = i 局 0 | 

i100- WT--. 2 避 妃 电 
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1980 1985 “1990 ,1995~2000 2005 2010 2015 1985 1990 1995 2000 2005 2010 2015 
年 份 Year 


图 5 南京 市 LAI 年 际 变化 趋势 
Fig.5 The interannual variation trend of LAI in Nanjing City 


如 图 5a 所 示 ; 从 全 市 整体 来 看 ,1988 一 2013 年 期 间 ,LAI<2 的 低 植被 覆盖 区 域 不 断 增 加 ,LAI>3 的 高 植被 
履 盖 区 域 面积 有 先 减 少 后 增加 的 趋势 。LAI<2 的 区 域 面 积 由 335.81 km 升 至 798.87 km ;LAI>3 的 区 域 面积 
由 203.97 km 减少 至 71.41 km (1996 年 ) ,而 后 又 增加 至 270.15 km ,与 1988 年 基本 持平 。 从 20 世纪 80 年 
代 开 始 ,全 市 城区 面积 不 断 扩张 ,土地 利用 动态 变化 数据 显示 从 1986 年 到 2013 年 间 , 居 民 点 及 建设 用 地 面积 
从 11.85% 上 升 为 27.2% ,耕地 及 林地 面积 不 断 减少 ,从 78.19% 下 降 到 60.84%。 低 植被 覆盖 区 面积 不 断 增 加 ， 
植被 遭 到 严重 破坏 ,高 植被 覆盖 区 域 的 面积 迅速 减少 ;进入 21 世纪 后 ,城市 面积 还 在 不 断 扩 张 , 低 植被 覆盖 区 
面积 继续 增加 ,而 由 于 国家 政策 以 及 人 们 意识 的 改变 ,在 加 速 城市 建设 的 同时 注重 保护 和 恢复 植被 ,高 植 
被 覆 善 区 面积 不 断 增 加 ,但 是 耕地 面积 不 断 减 少 。 

李 明 诗 等 以 南京 市 主 城区 为 例 ,以 10 期 (1986 至 2011 年 ) Landsat TMZETM+ 遥 感 数据 为 信息 源 , 采 用 线 
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性 光谱 混合 模型 提取 植被 覆盖 度 ,对 研究 区 植被 覆盖 时 空 动态 模式 特征 进行 了 评价 。 为 了 对 比分 析 李 明 
诗 等 评价 结果 ,本 文选 取 南 京 市 主 城区 (玄武 区 .鼓楼 区 ,秦淮 区 . 建 邮 区 .雨花台 区 栖霞 区 ) 统 计 分 析 时 间 序 
列 (1988 一 2013 年 )LAI 的 分 布 变 化 情况 。 图 5 表明 主 城区 1988 一 2013 年 间 ,LAI<2 的 区 域 不 断 增 加 ,LAI>3 
的 区 域 不 断 减 少 。LAI>3 的 面积 由 50.49% 降 至 35.75% ;LAI<2 的 面积 由 27.58% 升 至 42.20% 。 其 变化 趋势 
与 李 明 诗 ”等 得 到 的 南京 城市 绿地 时 空 动态 变化 一 致 。 其 主要 原因 是 1988 一 2013 年 南京 市 主 城区 城市 化 
快速 发 展 , 城 市 人 口 迅 速 增加 ,目前 全 南京 市 常住 人 口 总 量 为 818.78 万 人 ,常住 人 口 城镇 化 率 达 到 了 80.5% ， 
大 力度 的 房地产 开发 ,城镇 用 地 不 断 扩张 ,良好 的 植被 遭 到 破坏 ,导致 LAI>3 的 面积 下 降 , 而 随 着 生态 文明 的 
建设 , 林 草 植被 恢复 率 又 得 以 改善 ,导致 LAI 又 表现 出 一 定 的 波动 性 。 


3 结论 与 讨论 


本 文 应 用 BP 神经 网 络 模型 对 南京 市 1988 一 2013 年 时 间 序 列 的 叶 面 积 指数 进行 反 演 ,分 析 BP 神经 网 络 
模型 反 演 时 间 序 列 LAI 的 精度 情况 ,对 比分 析 LAI 的 年 际 变 化 与 土地 覆盖 变化 的 关系 。 得 出 如 下 结论 . 

(1) 模 型 拟 合 度 较 高 ,结合 2009 年 和 2010 年 的 实测 数据 ,两 年 的 平均 相对 误差 均 方 根 误差 .相关 系数 分 
别 是 0.2395 和 0.2174 .0.2962 和 0.2581 .0.7713 和 0.6844 ,各 项 精度 评价 指标 均 较 好 。 

(2)1988 一 2013 年 期 间 ,南京 市 LAI<2 的 主要 分 布 在 建筑 物 比 较 密 集 植被 稀 玻 的 城区 ,包括 江宁 区 北 
部 .六合 区 南部 浦口 区 东南 部 .栖霞 区 以 及 市 辖区 ;LAI>4 的 主要 分 布 在 紫金 山 . 老 山 .幕府 山 ,将军 山 . 铜 山 
等 植被 比较 密集 的 山区 。 

(3) 统 计 分 析 去 除 耕 地 后 南京 市 LAI ,在 1988 一 2013 年 期 间 , 低 植被 覆盖 (LAI<2 ) 面积 不 断 增加 ,高 植被 
覆盖 区 (LAI>3) 面积 先 减少 后 增加 ,耕地 面积 不 断 减 少 ,符合 南京 市 的 发 展 变化 规律 。 

(4) 统 计 分 析 南 京 市 主 城区 LAI,LAI 的 年 际 变化 与 其 他 学 者 得 到 的 南京 市 植被 善 度 变化 趋势 一 致 , 反 演 
结果 的 时 序 性 较 高 。 研 究 结果 表明 ,BP 神经 网 络 模型 对 时 间 序 列 LAI 反 演 具有 较 高 的 精度 ,为 实现 土壤 侵蚀 
模型 中 植被 覆盖 管理 措施 因子 C 的 定量 遥感 反 演 及 区 域 土壤 侵蚀 定量 遥感 监测 提供 了 新 途径 。 

但 是 ,由 于 BP 神经 网 络 模型 的 训练 存在 一 定 程度 的 误差 ,而 且 反 演 区 域 的 面积 较 大 ,植被 类 型 与 植被 群 
落 结 构 复 杂 多 样 ,农作物 的 种 植 状 况 复 杂 ;导致 时 间 序 列 LAI 还 感 反 演 精度 有 待 进一步 提高 , 反 演 方法 还 需要 
进行 进一步 研究 。LAI 高 度 概 括 植被 水 平 履 盖 与 垂直 结构 特性 , 利用 多 角度 遥感 数据 及 物理 模型 ,建立 多 角 
度 LAI 的 反 演 方法 ,构建 不 同 植被 结构 下 LAI 与 土壤 侵蚀 的 耦合 模型 ,以 及 多 源 遥 感 影像 定量 融合 ,各 种 LAI 
反 演 方法 交叉 使 用 ,为 区 域 土壤 侵蚀 定量 遥感 监测 提供 新 途径 ,将 是 未 来 研究 趋势 。 
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